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INTERPOLATION INEQUALITIES IN PATTERN FORMATION 


ELEONORA CINTI AND FELIX OTTO 

Abstract. We prove some interpolation inequalities which arise in the analysis of 
pattern formation in physics. They are the strong version of some already known 
estimates in weak form that are used to give a lower bound of the energy in many 
contexts (coarsening and branching in micromagnetics and superconductors). The 
main ingredient in the proof of our inequalities is a geometric construction which 
was first used by Choksi, Conti, Kohn, and one of the authors in [3] in the study of 
branching in superconductors. 


1. Introduction 


In this paper we establish some interpolation inequalities which are connected with 
the study of certain physical phenomena. More precisely, the inequalities that we 
prove are the strong versions of some already known interpolation estimates in weak 
form that play a crucial role in the study of pattern formation in physics. 

In many physical phenomena described by a variational model, in order to un¬ 
derstand why certain patterns are observed, it is natural to study the features of 
pattern with close to minimal energy. This requires a good understanding of at least 
the scaling of the minimal energy in terms of the model parameters. Upper bounds 
on the minimal energy are obtained by physically motivated trial patterns (Ansatz). 
Typically, Ansatz-free lower bounds rely on suitable interpolation inequalities that 
involve some functional norms related to the energy under consideration. 

The first two interpolation inequalities that we present here involve the BV-norm 
and the R'^-norm of a function u. The first estimate holds in any dimension d and it 
is stated in Proposition 0 below. The second inequality (see Proposition |1.2[ ) holds 
in dimension 2 for functions bounded below, and it improves the result in Proposition 


1.1| by a logarithmic factor. 

Here with H~^ we refer to the homogeneous negative Sobolev space H~^. Let u be 
a periodic function defined on the torus [0, A]'^ and with vanishing average. We recall 
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that the H ^-norm of u is dehned as follows 


tu 


H 


-1 = lliv 


-M\l 


:= inf 

j 


'[0,A]" 


|j|^ I j is periodic and V ■ j = u 




|V(pp where — A(p = u. 


It can also be dehned via Fourier transform: 


ll|V|-^«||^ = j m-"\nn)\Ydk, 

where / dk has to be interpreted in a discrete sense. 

Proposition 1.1. Let u : [0, A]*^ M. be a periodic function such that = 0. 

Then, there exists a constant C > 0 only depending on d such that 


\u\\4 < C'||Vn||f|||V 




( 1 . 1 ) 


Let us comment on the nature of (1.1), which is the model of three further in¬ 
terpolation estimates appearing in this paper. Its interest for applications is that it 
relates: 

• the L^-norm of Vm, with its obvious geometric interpretation via the coarea 
formula and its ensuing connection to an interfacial energy; 

• the L^-norm of with —Aip = u: J is a prototype for a held energy, 

e.g. the electrostatic energy of u interpreted as a charge distribution; 

• some L^-norm of u itself, which typically can be estimated below using the 
non convex features of the model where u plays the role of an order parameter. 


Let us go one step deeper: (1.1) is one end-point estimate in the one-parameter 
family of Gagliardo-Nirenberg estimates 

||m|Ip < C'||Vm||||||V|“^m|||, l<g<oo, (1.2) 

where the integrability exponent p on the left-hand side is determined by the relation 

1 _ 1 1 11 
p “ 2 ' g ^ 2 ' 2’ 

that is p = This exponent is the only one compatible with rescaling of length. 
At the same time, the estimates (1.2) are “extensive”, that is, compatible with taking 
volume averages, i.e. passing from dx to A“'^ dx, as can be seen from 


rewriting (1.2) as 


'[0,A]“ 


ujPdx < C 


'[0,A]‘* 


|V«|' 


1 1 
2'? 


'[0,A]“ 


1-1 


U\ 


1 1 
2 ■ 2 
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This extensivity is a crucial feature in the applications: the estimate should be com¬ 
patible with taking volume averages, since it should be oblivious to the artihcially 


introduced period A. In this sense, the family of estimates (1.2) is “orthogonal” to 
Sobolev estimates, which are saturated by localized functions. 


Among the family of estimates (1.2), the one with q = 2 (and thus p = 2) is 


particularly simple and can be established by Fourier representation. The estimates 
for 1 < p < oo are fairly easy to prove using Calderon-Zygmund theory, i.e. the 
maximal regularity for —A in L'^-spaces. The latter is no longer available for the 
end-point q = 1 and thus it is not surprising that this estimate was only recently 
established. It was hrst proved by Cohen, Dahmen, Daubechies, and DeVore (see 
Theorem 1.5 in i) using the wavelet analysis of the space BV, while our proof uses 
a technique introduced by Ledoux (see Theorem 1 in [I3]) to give a direct proof of 
some improved Sobolev inequalities. 


The following proposition improves the result in Proposition (1.1) by a factor In^ u 
and it holds in dimension d = 2 for functions bounded below. 

Proposition 1.2. Let u : [0, A]^ —>■ M 6e a periodic function such that u > —1 and 
/[o,A]‘* “ = 0- 

Then, there exists a constant C > 0 such that 


lluln^ max{M,e}||4 < C||Vm||i |||V| (1-3) 

The main ingredient in the proof of Proposition |1.2|is a geometric construction (see 


Lemma 4.2), which was first used by Choksi, Conti, Kohn, and one of the authors in 


[1] in the context of branched patterns in superconductors. 

These two interpolation inequalities are naturally connected to many physical prob¬ 
lems (coarsening, domain branching in ferromagnets, superconductors, twin branch¬ 
ing in shape memory alloys) whose energy is given by the competition of two main 
terms: an interfacial energy (described by a BV-norm) and a held energy (described 
by an iP“^-norm ). In Section 2 we will explain how estimates (1.1) and (1.3) are 


used to give lower bound for the energy in two different problems: coarsening and 
micromagnetics. 

Our last two interpolation inequalities are connected with the study of branching 
in superconductors and they both involve the Wasserstein distance. 

We recall (see for instance [20]) that the Wasserstein distance between two measures 
with densities u and v is given by 


W^iu, v) := inf 



\x - y\‘^dn{x,y)\ j d7i{-,y)=u, j dn{x,-)=v 
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(1.4) 










The measure on the product space vr is called transportation plan and it is admis¬ 
sible if its projections to hrst and second coordinates are measures with densities u 
and V respectively. 

A useful property of the Wasserstein distance is the Kantorovich duality (see Chap¬ 
ter 5 in 1201 ), which allows to write the Wasserstein distance in a dual form in the 
following way: 


v) = sup < / u{x)(p{x)dx + / v{y)il){y)dy \ >p{x) + ip{y) < |x — yY 


(1.5) 


This property will be useful in the proof of Propositions 1.3 and 1.4 below. 

Loosely speaking, the Wasserstein distance between u and n is a nonlinear version 
of a negative norm oi u — v. The hrst proposition states an interpolation estimate 
with the iL“^-norm of u (or rather of n — 1) replaced by W 2 {u, 1). 


Proposition 1.3. Let u : [0,A]'^ —)■ M 6e a periodic function such that u > 0 and 
Then, there exists a constant C > 0 only depending on d such that 

2d ^ d 

IKn-Cj + ll^+M < C'(llVnlli)^ (W|(n,l))5+3a, (1.6) 

3d 

where {u — C)+ := max{-u — C, 0}. 


In the recent paper [H], M. Ledoux showed that our inequality (1.6) is actually a 


particular case in a all family of interpolation inequalities (called Sobolev-Kantorovich 
inequalities) that hold in the more general setting of non-negatively curved (weighted) 
Riemannian manifolds. The proofs in [H] rely on the use of heat hows and Harnack 
inequalities. 


Let us comment a bit on interpolation estimate (1.6). Because of the formal scalings 


||Vn||i has unit of length ^ x volume, 
hP|(n, 1) has unit of length^ x volume. 


{u - C). 


has unit of 


volume p, 


the pair of exponents on the right-hand side of (1.6) is determined by the integrability 


exponent p = appearing on the left-hand side. The pair of exponents is the only 
one compatible with rescaling of length and taking the volume average. The exponent 
p could not be inferred by a simple scaling argument. It is the exponent that appears 
on the left-hand side in the linear interpolation estimate 


|'U|| 2+3d < CllVtill^ 


2d 

2+3d I 


2d 


d+2 

d 
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In this sense, 1), at least as seen from the “peaks” of u (i.e. the part of u that 

is much larger than 1), behaves as the negative fractional Sobolev norm 


2d 


d+2 

I d _ 

I d+2 — 

d 


2d d+2 

V d . 


The last nonlinear interpolation estimate of this paper replaces the negative norm 


in Proposition 1.1 and the Wasserstein distance in Proposition 1.2 by a mixture of 


both. We recall that by the iP-method of interpolation spaces (see Chapter 3 in 0 ) 
a norm || ■ || intermediate between some norms || ■ ||(o) and || ■ ||(i) can be obtained via 
the construction 


|m|| = sup inf 


u - n||(o) + /i » ||n 


By the equivalence of approximation and interpolation theory (see Chapter 7 in 0 ) 
in term of 

||m|| ~ ( sup inf {||m — n||(o) I ||n||(i) < ) , 

\M>0 / 

we also obtain a representation with mixed homogeneity 


\u\ 


sup inf < ix\\u - v\\fo) + (jx <> ||n||(i) 


2 - 6 > 

2 


In Proposition |1.4[ with an application described in the next section in mind, we 
replace ||m —n||(o) by IV^iu, v) and ||n||(i) by the homogeneous fractional Sobolev norm 


|V| 2 n|| 2 , and choose 9 = This leads to 

sup inf < fiW 2 {u,v) + 

^>0 1 V / 


d+1 
2 'i d+3 


(1.7) 


= sup inf \i'd+iW 2 {u,v) + i' ‘*+i|||V| 2 n ||2 

i/>0 


d+1 

1 < I 2 1 d-\-3 


Proposition 1.4. Let u : 


-)■ 


he such that u > 0 and J u < oo. 


Then there exists a constant C > 0 only depending on d such that 


|m||m± 3 < CII■ sup inf \i"i+^W 2 {u,v) + u ‘*+i|||V| 2 n ||2 

3d+i u>0 ^ 


1 

2 I 3 


In our proof, we will use the following representation of the fractional Sobolev 


norms 






Let us comment on Proposition We claim that it can be interpreted in the 
sense of acts as |||V|“^(-u — u)||i”. Indeed, if in (1.7) we replace W^iu^v) 

by II|V|“^(m — u)||i, we have 


sup inf </i|||V| ^(M-u)||i+fp ^^|||V| 2 ^ 112 ) 

^>0 1 V / 


d+1 
2 ^ d+3 


According to the above discussion, this expression is equivalent to 


sup inf 

^>0 


-2 


(u-u)||i+/i-^|||V|-W|h , 


where 9 = The norm of this interpolation space [H 


-2 




0 has the same scaling 


as the norm of the fractional Sobolev space H 


(i_e){_2)+0(-i) _ .- 


(i_e)i+ei 


2d 

_ TJ d+3 

— d+3 ! 

d+1 


that is 


V 


2d 

d+3 11 


d+3 . 
d+1 


With this substitution, the nonlinear interpolation estimate of Proposition |1.4 turns 


into the much more standard linear interpolation estimate 


u 


3d+3 

3d+l 


2d 

< C'||Vu||f+^|||V 


2d 

d+3 'll 


d+3 

3d+3 

d+3 

d+1 


In this sense, W 2 {u,v) acts as || |V|“^(m — u)||i. 

When deriving an Ansatz-free lower bound in the physical applications, only a weak 
version of these four interpolation inequalities is needed, in the sense that it is enough 
to have these estimates with the L^-norms on the left-hand side replaced by the weak 
L^-norm. The contribution of this paper is to pass from the weak formulation to the 
strong one. 

The paper is organized as follows: 


In Section 2 we expose the connection between the above interpolation in¬ 
equalities and the study of pattern formation by considering three examples: 
coarsening, micromagnetics, and superconductors. 

In Section 3 we prove Proposition using a method introduced by Ledoux 
in [13]. 


In Section 4 we give the proof of Proposition |1.2[ whose main ingredient is 
the geometric construction given in Lemma 4.2, introduced in |1]. 


In Section 5 we prove Propositions 1.3 and 1.4 
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2. Three models: coarsening, micromagnetics, and superconductors 


In this section we expose the connection between our interpolation inequalities and 
three different physical phenomena. 

The first model that we consider describes thermodynamically driven demixing, 
the second one describes the magnetization in a ferromagnet. In the first model, we 
are interested in the phenomenum of coarsening of the spatial phase distribution, in 
the second model, we are interested in the phenomenum of domain branching in a 
strongly uniaxial ferromagnet towards the sample surface. 

In both of them a crucial ingredient in the proof of a lower bound for the energy is 
a weak version of our first two interpolation inequalities (Propositions 0 and |1.2[ ). 

The last model that we consider describes type-I superconductors. Also here, the 
phenomenum of interest is branching: it is the magnetic flux carrying normal phase 
that branches when approaching the sample surface. In this phenomenum, several 
different regimes occurs and Propositions 1.3 and 1.4 are related to the analysis of 
two of them. 

The discussion in this section is not rigorous; the main aim is to show how the 
quantities involved in our interpolation estimates arise naturally from the energy 
functionals associated to these three models. 

1. Coarsening. The Cahn-Hilliard equation models the thermodynamically driven 
demixing of a two-component system. The relative concentration is described by an 
order parameter, say on a large torus, u : [0,A]'^ —)■ M. The driving free energy is of 
Ginzburg-Landau form 


m = 2 


[0,A]“ 


+ {i-uy)dx. 


( 2 . 1 ) 


The postulate that the free energy decreases while the order parameter is conserved 
in time (i.e. ^ /[q a]'* = 0) is satisfied by the evolution equation 


dS 

dtu — A—— = dtu — A(Am + u — u^) = 0, 
ou 


( 2 , 2 ) 


which is the Cahn-Hilliard equation. 

Starting from a small perturbation of the critical mixture u = 0, numerical simu¬ 
lations show after an initial stage the emergency of two convoluted regions in which 
u is very close to either its equilibrium value 1 and —1, respectively. These regions 
are separated by a shortly modulated characteristic transition layer given by the one¬ 
dimensional minimizer of (2.1). In the sequel, the evolution is essentially geometrical 
(named after Mullin and Sekerka) and is driven by a reduction of the interfacial 
energy. This leads to a coarsening of the two complementary domains, i.e. to an 






increase of their average length (as embodied by its average width or radius of curva¬ 
ture). Eventually, this coarsening stops because the average length-scale reaches the 
artihcial period A. The coarsening is observed to be “statistically” self-similar in the 
sense that e.g. the two-point correlation function is (approximately) self-similar. It 
is characterized by the exponent 1/3 in the sense that average length scale as 

In |T2], Kohn and one of the authors proved quantitatively that, independent of the 
initial data provided that they are well-mixed, coarsening cannot proceed at a faster 
rate. The strategy for the proof has since been applied to many models that feature 
coarsening (see e.g. the discussion in [3]) and relies on the gradient flow structure. 
Morally speaking, it consists in converting information on the steepness of the energy 
landscape (how quickly does the energy decrease as a function of the intrinsic distance 
to the well-mixed state) into information on the speed of relaxation (how fast does 
the energy decrease as a function of time), see the discussion in |T7]: if the landscape 
is not too steep, the relaxation is not too fast. 


As can be easily inferred from (2.2) the Cahn-Hilliard evolution is the gradient flow 


of the Ginzburg-Landau energy £ with respect to the Euclidean geometry given by 
the inner product in H~^. Hence showing that the energy landscape is not too steep 
means establishing an estimate of the form 

1 


T(u)|||V 


-1 


m ||2 > 


(2,3) 


the power of the system volume A'^ being dictated by the need to take volume averages. 

We now heuristically point out the connection to the interpolation estimate from 
Proposition 1.1 As mentioned above, the evolution relevant for coarsening is well 


approximated by the free-boundary problem named after Mullin and Sekerka, which 
is the gradient flow of 


^Msiu) — 


CoJ\Vu\ ifue{-l,l} 


-l-CXD 


else 


with respect to H Hence replacing £ by £ms in (2.3) leads to 


llVulh 


IV7I-1 II ^ A I 1 11 ,, 

|V| «l|2>pA. = -||t.||., 


which is the squared version of our interpolation estimate. 

Let us now turn to Proposition 1.2[ If the mixture is strongly off-critical, by which 
one understands that 








the above strategy gives the correct scaling of the prefactor in the upper coarsening 
bound in the volume fraction (j) in dimensions d > 2. In low volume fraction, numerical 
simulations show that the region covered by the minority phase u = 1 breaks up into 
many connected components which quickly relax to have approximately round shape. 
Coarsening then proceed by Ostwald ripening: the large balls grow at the expense of 
the small ones that eventually collapse. 

Asymptotic analysis for d = 2 shows that there is a logarithmic correction to the 
prefactor in the coarsening rate [16]. It is the interpolation estimate in Proposition 


1.2 that correctly captures this logarithmic correction, as shown in [5|, where also the 


weak form of (1.3) was established 


2. Branching in Ferromagnets. We are interested in a ferromagnet in form of a 
slab of thickness 2t] for simplicity, we assume periodic boundary conditions in the 
horizontal directions so that the fundamental domain is given by [0, A]^ x {—t, t). The 
magnetization m is a unit-length vector held m : [0,A]^ x —)■ that locally 

indicates the mean direction of the spin. 

The ground state minimizes the Landau-Lifshitz energy functional (as modihed by 
Brown): 

E{m) = d^ / \Vm\^dx+Q / {ml+ml)dx+ \\V\^^'V-m\'^dx. 

J[0,A?x{-t,t) J[0,A?x{-t,t) J[0,A?xR 

(2.4) 

The hrst term, the exchange energy, is of quantuum mechanical origin, and models 
a short range attraction between the spins and comes with its intrinsic length scale 
d. The second term, the anisotropy energy, comes from the interaction of the mag¬ 
netization with the lattice structure of the metal; here, it is uniaxial and favors the 
third axes and thus the values m = ±(0, 0,1). The last term, the stray held energy, 
can also be written directly in terms of the stray held h, which is determined by the 
static Maxwell equations: 


'[0,A]2xI 


\h\‘^dx, where V ■ (h -|- m) = 0, V x h = 0. 


(2.5) 


Both equations are to be interpreted distributionally in where m is extended 
trivially outside the slab. In other words, h is the Helmoltz projection of the extended 
m. There is a mathematical analogy from magnetostatics convenient for the intuition: 
there are two sources for the stray held, namely ’’volume charges ” coming from —V -m 
in the bulk x {—t,t) of the sample and ’’surface charges” ± 777.3 at the surface 
X {—t,t} of the sample. Like in electrostatics, both give rise to the potential held 
h. 
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Following Hubert (see Chapter 3 in m), let us euristically explain why branching 
occurs; see also [IHIIIH]- A constant magnetization in direction of the easy axes, say 


m = (0, 0,1), would have zero exchange and anisotropy energy. However, by (2.5), we 
would have h = m so that the stray held energy per area of cross section would be 2t. 
However, the if“^-norm of distributional divergence V -m, which here comes from the 
ma-component at the sample surfaces [0,A]^ x (that is the surface charges), 

can be reduced by horizontally alternating between m = ( 0 , 0 , 1 ) and m = ( 0 , 0 , — 1 ), 
and thus alternating the sign of the charges. In fact, for magnetization that only 
depends on x' = {xi^X 2 ) and alternates with a period u <^t, the stray held energy 
would behave as the squared H~^-noTm in the horizontal variables x' = {xi,X 2 ) of the 
surface charges m 3 , that is like \ \V\~^m 2 \'^dx'. The exchange energy prevents 

such jumps of m, but they can be replaced by smooth xi-dependent transition layers 
(Bloch walls) between the magnetizations (0,0,—1) and (0,0,1). These transition 
layers are determined by a balance between exchange and anisotropy energy in form 
of 




dm 


dxi 


dxi + Q {ml + ml)dxi 


= d^ 


1 — m3 


dm 3 

dxi 


dxi + Q / (1 — m\)dxi 




dm3 

——dxi = aq) 2. 
dx\ 


Dimensional analysis shows that the width of the transition layers scales as 
^waii ~ dQ~^ and the above inequality shows that the wall energy per area in vertical 
direction is given by dQ^. This suggests to consider the reduced energy functional 


E{m) = dQ^^ / \Vm 3 \dx+ / ||V| ^c^smsl^dx (2.6) 

i[0,A]2x(-i,p i[0,A]2xM 

subject to m 3 = ±1. From this reduced functional it is clear that there is an optimal 
period codomain for up-down domains alternating in the a:i-direction and separated by 
Bloch walls. The wall energy of such a conhguration per area in the cross section X 1 -X 2 
behaves as {specific wall energy) x {wall area per cross section) = dQ^ x 
From the above remark on the if“5-norm we have by a dimensional argument that 
for oodomain ^ the stray-£eld energy scales as Udomain- Hence the optimal period 
scales as uidomain ~ {dQH)^, which is consistent with Uwaii ^ (^domain -C t provided 
dQ^ t. The energy per cross-sectional area scales as {dQH)^, beating the uniform 
magnetization in the above regime. However this is not the scaling of the minimal 

energy! Intuitively, it is advantageous to have the period UJdomain depending on 0 : 3 : 
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for ~ =bt, the domain width uJdomain{,x^) should be vanishing in order to have mean 
cancellation in the ’’surface charges” given by m^\ for 0:3 away from the surfaces, 


the domain width uidomainixs) should be large in order to avoid wall energy. Such a 
height-dependent domains width can be realized by domain branching. However, 
there is no free brunch: the branching of domains means that the interfacial layers 
between ( 0 , 0 , 1 ) and ( 0 , 0 ,— 1 ) are no longer vertical but tilted, hence they carry a 
(volume) charge V ■ m and thus generate a stray held. If this tilt is small, at least in 
the bulk, this suggests to consider the following anisotropic reduction of ( 2 . 6 ): 


E{m) = dQ^ 


'[0,A]2x(-t,t) 


iV'm^ldx 


'[0,A]2x]R 


V'l ^dsmsfdx, 


(2.7) 


subject to m 3 = ±1. Note that with our convention m 3 = 0 for |a; 3 | > t, the hniteness 
of the anisotropic version (2.7) of the stray held energy implies that m 3 ^ 0 for 0:3 —>■ 
±f, that is, inhnite branching. Now thanks to its anisotropic character, the reduction 
(2.7) has a scale invariance that allows for the following non-dimensionalization 


x' = {dQ 2 ) 3 t 3 x' (and thus A = (dQa) 3 ^ 3 A), 

X 3 = tx 3 , and -^E = {dQ^)H^ ^E, 
A^ A2 


( 2 . 8 ) 


that removes all parameters. This suggests that the domain width in the bulk scales on 
average as uidomain ~ {dQ^)H^, and the energy per cross-sectional area as {dQ^)H^, 
which indeed beats the unbranched case. 


The passage from (2.4) to (2.7) has been made rigorous in [TH] on the level of a 
T-convergence result for inhnite cross-sectional area. A key ingredient for this is a 
lower bound on the minimal energy per cross-sectional area, i.e. ^mini7, that is 
independent of the artihcial periodicity A S> 00 domain- It is here that the interpolation 
estimate from Proposition 1.1 comes in. Let us show how, w.l.o.g. on the level of 


the non-dimensionalization (2.8). By Poincare’s inequality in X 3 (recall that m 3 is 
supported in (—1,1)), Young’s inequality and the interpolation inequality in x' (and 
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thus for d = 2), we have as desired: 


E{m) = 


'[0,A]2x(-l,l) 
rl 


iV'm^ldx + / ||V'| ^d^m^l'^dx 


'[0,A]2xR 


> 

rsj 


> 

rs_/ 


> 


iV'rrisldx' + / ||V| dhs 


'-1 VJ [ o , a ]2 


'[0,A]2 


'-1 V^[0,A]2 
rl r 


\V'm'i\dx' 


VI ) dxs 


'[o,A]2 


J-l ~'[0,A]2 

> 


Imsl sdh'dxa 


3. Branching in Superconductors. 

The so-called intermediate state of a type-I superconductor is characterized by 
penetration of the magnetic held in some parts of the material which leads to the 
formation of normal and superconductive domains. The superconductive regions are 
characterized by the expulsion of the magnetic held, this is the so called Meissner 
effect. In mu the mathematically rigorous study of this pattern formation prob¬ 
lem, via energy minimization, was initiated. In particular, the authors established 
rigourous upper and lower bounds for a reduced energy associated to this problem, in 
diherent regimes depending on the value of the external magnetic held (small, inter¬ 
mediate, or close to critical). Similarly to the case of micromagnetics, interpolation 
inequalities can be seen as playing a crucial role to prove an ansatz-free lower bound. 
In the sequel we briehy describe the physical model. 

Following [9], for A,f > 0, we consider a wave function -0 : [0,A]^ x {—t,t) —)■ C, 
which plays the role of the order parameter, and a vector potential A : [0 ,A]^ x 
{—t,t) —t We introduce the following quantities: 

the density of superconducting electrons 


the magnetic held B := V x A, 

the covariant derivative VA'f '■= VV' — iAf), 

the kinetic energy iVyiV’T) 

and the superconductive current j := X{'ip'VA'f)- 

We observe that, if is written in polar coordinates if = pe*®, then 

= |Vp|^ + p^|V0-A|2 and ] = - A). 
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Sometimes we will use the notation B' to denote the hrst two components of the 
magnetic held: B = {Bi, B 2 , B^) = {B',B^). 

There are three parameters that govern the behavior of the material: the external 
magnetic held -Bexo the coherence length which measures the typical length on 
which -0 varies, and the penetration length A, which describes the typical length 
on which the magnetic held penetrates the superconducting region. The Ginzburg- 
Landau parameter is given by k = |. 

For any pair such that physically observable quantities p, B, j are [0,A]^- 

periodic, we dehne the Ginzburg-Landau functional: 


Eo{'iP,A) := [ [ \B - B,,t\^d(^.9) 

J[0,Apx(-t,t) \ ^ J ./[0,A]2xR 


[0,A]2x(—t,i) 


K 


'[0,A]^x]R 


IVpr + P I+ -(1 - py dx 


\B - Bext?dx, 


( 2 . 10 ) 


where i?ea;t = (0, 0, <F) is the external magnetic held and the penetration length A is 
normalized to be 1 (that is, in this unit, k represents the inverse of the coherence 
length). The Meissner ehect, namely the fact that the kinetic energy disfavors the 
magnetic held where the material is superconducting (i.e. p > 0) can easily be 
seen from this formula: In simply connected regions where p^ is positive, the vector 
potential A wants to be close to gradient (the gradient of the phase 6) so that the 
magnetic held B wants to be small. 

The type-I superconductors correspond to the regime of k small {k < \/2). In this 
regime, there is a positive surface tension that leads to the formation of normal and 
superconductive regions corresponding to p ~ 0 and p ~ 1, respectively, separated 


by interfaces. Indeed, using an identity on fhe energy Eq in (2.9) can be 

bounded below by (see [9], Lemma 2.3) 


Eo{'ijj,A) > 


'[0,A]2x(-t,p 


K 


V2 


1-^ |Vp|^+ i?3-^(l-/) 


K 


K 


B'Y - + —$ 

72 


dx + 


72' 

B - Bext^dx. (2.11) 


'[0,A]2x(R\(-t,t)) 
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We observe that, under the sharp Meissner condition p^B = 0, the sum of the first 
two terms on the right-hand side can be written as 




K \ l2 / ^ « 




dx 


'[0,A]2x(-i,p 


(^1 “ + yX{p>0}(l - + X{p=0} ^ 


This is a Modica-Mortola type functional with a degenerate double-well potential 
given by 

^(p) = X{p>o}(l 

After introducing the new order parameter y = 1 — and rescaling according to 
X = ^x, B = -^B, $ = a Modica-Mortola type argument leads to the reduced 
functional in the regime k 1: 




-\Wx\ + \B'\^ + iBs-xr]dx + 


'[0,A]2x(R\(-t,p) 


l-B — (0,0, ^)\'^dx 


( 2 . 12 ) 

subject to the constraints: y G {0,1} (y = 0 corresponds to the superconducting 
phase, and y = 1 to the normal phase). 


V ■ B = 0 everywhere. 


and 

(1-X)5 = 0 in [0,A]^ X (-f,t). 


where V ■ B = 0 comes from the Maxwell equations in form of 5 = V x A, while the 
last constraint {1 — x)B = 0 represents the Meissner effect. 

2 

We now perform the following anisotropic rescaling: x' = tsx', x^ = tx^, B' = 
t~^B' inside the sample [0, A]^ x (—t,t), and x = t^x (and thus A = f^A), B = t~^B 
(and thus <h = t“5<l>) outside the sample. If we dehne z/ = we get the energy 
functional in the regime t 3> 1 (for simplicity of notations we drop the " on functions 
and variables): 


E{x,B) = 


'[0,A]^x(-l,l) 


^\Vx\+x\BT]dx+- 


^ d[0,A]2x(R\(-l,l)) 


\B — (0,0, <h)p(ia;, 
(2.13) 
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dx. 



subject to the constraints 

X ^ {0,1}, 

+ V' ■ ixB') = 0 in sample, 

X = B 3 at surface, 

V • -B = 0 outside sample. 

Note that in this regime of t !:§> 1, the penalization of dsx fades away, while the 
penalization of B^ — x turns into the hard constraint i ?3 = y in the sample, so that 
together with the Meissner effect 5(1 — y) = 0, V ■ 5 = 0 turns into the transport 
equation 

dsx + V-ixB') = 0. 

As for the case of micromagnetics, also in this case a branched pattern is observed: 
X alternates between the two phases x = 0 and x = 1 011 a length-scale which 
decreases while approaching the boundaries { 0:3 = ± 1 }. 



Figure 1. 

Let us have now a closer look at the quantities involved in the energy E{x,B). 

We hrst observe that, thanks to the last two constraints in (2.14), the last term in 
the energy is estimated by 


'[0,A]2x(R\(-l,l)) 


\B — (0, 0, 

15 


> 


|V| - <*)||2, 


(2.15) 






























since (cf. ( 1 . 8 )) 


|||V|-2/||2 = 27rinf / \B\^dx : V ■ i? = 0, B,{;0) = f 

L d [0,A]2x (O,oo) 

Here, we denote by Ihe weak limit of y as X 3 t 1; in particular, may 

not be a characteristic function. 

We consider now the second term a] 2 x(_i 1 ) x\B'\'^dx: Using again that y satishes 
the continuity equation, cf. (2.14), 


dsX + V-ixB') = 0, 

the classical Benamou-Brenier result in optimal transport theory [1] implies that for 
any slice { 0:3 = z} for z G (— 1 , 1 ) 




x\B'fdx > 


(2.16) 


Combining together (2.15) and (2.16), we deduce that there exists a slice z G (—1,1) 
such that 


E{x,B)> I |V'x| |dx' + lU|(x| ,X| ) + -|||V| Hx\ 


'[0,A]2 


U 


{2:3=1} 


(2.17) 

The quantities on the right-hand side are the ones involved in our Proposition |1.4[ 


In [HI Hj rigourous upper and lower bounds for the energy (2.13) are established. 
The study of minimizing conhgurations reveals different regimes, depending on the 
parameter i/ 1 and on the value $ of the external held. 

Here we list three different regimes with the corresponding behavior of the energy: 

(1) For 17 1 and 1 — <F ^ 1 the minimal energy per area behaves like: 

min U ~ (1 — $) Ins (1 — 4)). 

In this regime, the minority phase is not connected, as shown in Figure 2 
below. The interpolation estimate used for this regime is again the one with 


the logarithmic gain, established in Proposition 1.2 


(2) For u 1 and ur ^ ^ 1 the minimal energy per area behaves like: 


A ^ min U ~ 4) 3 . 


In this regime, the minority phase is connected and we have uniform branching, 
as shown in Figure 3 below. In this case, if we further simplify the model 

letting 17 —)■ 0 , the last term |||V|“^(x|^^ — 4 ))il 2 turns into the constraint 
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Figure 2. 


and in order to bound from below the energy -E(x, B) it is enough 
to bound from below the quantity: 



Vx\dx' + W^{X,^)- 


Let us show how, for this regime, our Proposition |1.31 leads to the right scaling 
for the minimal energy. By (2.16) and Young’s inequality we have for some 
slice z G (—1,1): 


E{x,B) = 


qo,A]2x(-i,i) 


(\y'x\+x\Bf)dx> I + 


'[O.A]" 




> 


'[0,A]2 




(2.18) 


^l{r.3=z} 


We apply now our inequality (1.6) for d = 2 to the function u = — ^ j m 
order to get 


$ 


-1 


max{M —2, 0}|| 


\\Vu\\! W2(n,l)3 


\ II r r\ r\^ II q A 2^ — - 

> II — 9. nl-l ^ o 


r\j TV ^ 3 
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(2.19) 





































Combining (2.18) and (2.19), we obtain 



(3) For v <^1 and $ ^ z/f the minimal energy per area behaves like: 

hr^minE ~ <Fz/“7. 

This regime differs from the previous regime by the non-uniformity of 

(see Figure 4 below). In this last regime, after rescaling by a suitable power of 
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z/, the lower bound for the energy can be deduced from Proposition |2.1| below 


d = 2 reads (after taking power = !)• 


(applied in the case of dimension d = 2). More precisely, inequality (2.21) for 

3d+3 _ 

3(i+l 7^ 


'[0,A]- 


9 

{u — dx <C ^11 Vn||i + n) + || |V|“^(n — <h)|| 2 j . (2.20) 


We observe that if we set u = Mu, v = Mv (and thus <h = M<h), x = ix (and 


thus A = iA) the quantities on the right hand-side of (2.20) scale as follows: 


liVnIli =£M||Vh||i, 

Wi{u,v) = tMWi{u,v), 

|V|-5(n - ^)\\l = eM^\\\V\-Hv - <l>)|| 


We now choose i = M = u a, use the scalings above and multiply hy u 9 in 


(2.20) to deduce: 


u 7 / {u — i>)ldx < C l^\\Vu\\i + W2{u,v) + u ^|||V| 2(-0 —<h)|| 
'[o,A]2 


Using (2.17) and applying the above inequality with u = and v = 

A|{^ 3 ^i}) we get the desired lower bound: 


A ^minP^ > $z/ 7 . 


Proposition 2.1. Let u, v : [OjA]*^ —)■ M 6 e periodic functions with u, n > 0 and 
A“'^ f u = f V = ^. 

Then, there exists a constant C > 0 only depending on d such that for any z/ > 0 

, 3d+l 

with d) < f^iysd+s we have 


/ 3d+l \ 

yu — iy9d+3 j 


< C ("II Vm||i -|- z^'Z+i hU|(n, n) -1- zz-z+i II |V| 2 (y 

3d+3 \ 

3d-fl 


3d+l 

4)11^)’". 

( 2 . 21 ) 


We write here the interpolation estimates (2.21) in additive form, since this is the one 
useful for the application. We prefer to state Proposition IM in multiplicative form, 
since this is the standard form for interpolation inequalities. In Section 5 we will see 
that Proposition lA follows easily from Proposition |2.1 
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Figure 4. 


3. Interpolation inequality in general dimension 


In this section we give the proof of Proposition 1.1 We start by recalling the weak 


version of estimate (1.1). We remind the dehnition of the weak L^-norm of a function 
u: 

\\u\\n,-p := sup/r|{|M| > 

/ j .>0 
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Lemma 3.1 ([T^). There exists a constant C < oo only depending on d such that 
for all periodic functions u : [0, A]'^ —>■ M, with f u = 0, we have 


|m|L_| < C\\Vu\\i |||V 


l-l„ 


u\ 


This lemma is proven in [T^. Here, for the sake of completeness, we give the proof 
of this weak estimate, since it is also useful to prove the strong version (1.1). 


Proof of Lemma TJ_. For simplicity of notations, in the following we will write “a < 6” 
to mean that there exists a positive constant C only depending on d such that a < Cb, 
and f u to denote u{x)dx. By a scaling argument in x, it is enough to show 


supp5|{|M| > fi}\ < ||Vm||i + |||V| ^ 

U>0 

Indeed, the change of variables x = Lx yields 

sup/i5|{|M| > p}| < L"^||Vm||i + 

M>0 


uwl 


i-i„ 


u\ 


2 ) 


where the symbol V denotes the gradient with respect to the new variable x. The 


choice of L = HV^ili |||V 


-1 


u\ 


yields 


SUp/i3|{|M| > p}| < llVwllf 
M>0 


1-1 


u 


Raising to the power 3/4 we get, as desired. 




ivr^iif. 


(3.1) 


< iiv^iif 

For an arbitrary level /i > 0 we introduce the signed characteristic function 
of the /i-level set of u\ 

{ 1 for /i < M 

0 for -p < M < p 
— 1 for u < —/i 

We select now a smooth symmetric f){x) > 0 supported in {|a;| < 1} with J ipdx = 1 
and dehne the Dirac sequence fjR^x) = Consider then the mollihcation of a 

function v on scale R, dehned as vr := 'iJjr * v. We have the identity 

j Xf^U = j Ur) + j XwRU- 

Using the duality between and H~^ in the second term on the right-hand side, 
we get the inequality 


/i j \Xf^\ < j Xl^U< \\U - + II VX;,,r||2 
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l-l„ 


M 2- 


(3.2) 




On the one hand, since 'ipR is supported in {|x| < i?} we have 

\\u - UrWi < i?||VM||i. 

On the other hand, using the dehnition of XfM,R and of '0 _r, we deduce 


\\Vx,,Rh<\\^i^R\\i\\xj2 = R-W\\ 


\x^. 


Plugging (3.3) and (3.4) into (3.2), we get 


A- / |X„I <-R||Vu||i + fl-‘||Vi/.||i / lx 


ivr'ti 


2- 


The choice of i? = /i a thus yields after multiplication with /is; 


/ lx,.I < I|V«||, + ||V-V,||i / lx 


ivr'ul 


(3,3) 


(3.4) 


Using Young’s inequality, we may absorb the hrst factor of the second term on the 
right-hand side and obtain the desired estimate. 

□ 


We give now the proof of Proposition [T^ The interpolation estimate (1.1) was hrst 
established by Cohen, Dahmen, Daubechies, and Devore (see Theorem 1.5 in |6]) by 
wavelet methods. We give here an elementary proof, which uses an idea by Ledoux 


Proof of Proposition o By scaling in x as in the proof of Lemma |3.1| it is enough 
to prove 

/ M5<l|v«||i + |||V|-'«||l. 

For an arbitrary level /i > 0 we use the signed characteristic function Xfi dehned in 


(3.1). Following an idea of Ledoux [13] for the proof a similar interpolation inequality 
we introduce a large parameter M S> 1 to be adjusted later. We consider as before 
Xii,R = '0R * X^i^ where 'i/jr is the convolution kernel dehned in the proof of Lemma 
13.11 We have: 

j Xi^u = JiXfi- Xf^,R)u + j Xti,RU 

= / (X/x-XM,i?)“+ / {x^l - X^l,R)u + / 

J{\u\<M^} J{\u\>M^} J 
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Using that ||x/i — X^l,R\\oo < 2, we obtain the inequality 


/ 


< Mu j - Xi,,r\ + 2 

< M^R [ |Vx^| + 2 


/ l«l + / 

\u\ + / X/^,rU- 




We multiply with /x 3 and choose R = jj, 3 as in the proof of Lemma 
over fi G (0, cxd) and using the duality between and H^, we get 


3.1 


Integrating 


_ 2 
fi 3 


< M 


/*oo /» 


poo 

|Vx^|da:ci/r + 2 / /i"i 
Jo 

' _2 \ 




<M / I Vx^|(ia:(i/i + 2 / fi 3 


q|n|>M/i} 


h ^Xti,RdfJ- 


where we keep the abbreviation R = fi 3 . 
On the left-hand side we have 




u\dxdfi 


u\dxdfi 


_2 f f _2 /■ 4 

/ /i 3 / |M(a;)|(ix(i/i = / |M(a;)| / /i 3dfidx = 3 / |m|3. 

io i{|n(a:)|>/i} J Jo J 

We address the three terms on the right-hand side one by one. We start by the second 


one: 



/ 




^Ar“^|n(a;)| 


|M(a;)| 


fi ^dfidx = 3M 3 


'0 


4 

li 3. 


We now address the first term. By the coarea formula we get 



\'Vx^l\dxdfi = / (Per({M >/i})-|-Per({-u < —/i}))(i/i = II Vm||] 
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/i 3/i' 3 / yxn,R-^Xii',R'dxdii'dii 


where Per(y4) denotes the perimeter of A. Finally we consider the last term (with 

V I I 

rx3 /‘oo 2 

Jo Jo J 

= 2 f j f (-A)x^,r XRL'^R'dxd/i'd^ 

Jo Jo J 

= 2 / f [ i>R'*{-Ai>ji)*XiiXii’dxdiJ,'diJ, 

Jo Jo J 

/■“ F 2 _2 

<2/ / /r 3;x' llx^lli ||x^,||ood/i'd/r 

Jo Jo 

, foo i-fi ^ 2 

= 2||AV’||l / / /i“3/i' 3/^-2 

JO Jo 
roo /-M 2 

= 2||A?/>||i / / /i'"3dA i|X;,||id/i 

Jo Jo 


6||A^/’||i / = 6||AV’||i 

Jo 

^IIAV’lli f \u\t 



|li(x)| 


fxsdfxdx 


These inequalities combine to 


4 

m 3 


< M||Vm||i + 6M 3 


H5 + I^1|A.A||i [Ms 


ivr^niu. 


We obtain the desired estimate by absorbing the middle right-hand side term for M 
large enough and absorbing the hrst factor of the last right-hand side term by Young’s 
inequality. □ 


4. Proof of Proposition 11.21 


In this section we prove Proposition |1.2[ We begin by recalling a geometric version 
of estimate (1.3), which was established by Conti, Niethammer, and one of the authors 
in M. 


Lemma 4.1 ([H])- Let x ■ [0;^]^ {OA} be a periodic characteristic function with 
volume fraction <F := J y <C 1. 
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Then, there exists a constant C > 0 such that 


$ Int - < C ( A-2 
$ 


iVxlJ J • (4.1) 

The proof of this Lemma made use of the following geometric construction, that 


plays a crucial role also in the proof of our strong estimate (1.3). 


Lemma 4.2 (|S])- For any periodic function x ■ [0, A]^ —)■ {0,1} and R L there 
exists a potential ((>r,l{x) G [0,1] such that 


X 


< 


R / I^X'I + / X't'R.L, 


j max{-A.^R,L,0} < ^ ' 1 j J x, 

f 't'R.L < L^R-^ f X- 


(4.2) 

(4.3) 

(4.4) 


We note that for L = i? we could just choose = 'ipR *x = XrI the interest here 


is the logarithmic gain In ^ A for L ^ R 


R 


Remark 4.3. We observe, for later reference, that for any function (f)'{x) G [0,1] we 
have 

uA' 

R- 

Indeed, we have 


j ■ V0' < R-^ (^In-i j 


X- 


(4.5) 


V (pR^L ■ V <p' 


= J {-A(j)R^L)(p'< J max{-A0R,L,O}0' 

< j max{-A(j)R^L,0} < R~‘^ ^ J x, 

where in the hrst two inequalities we have used 4>' > 0 and 0' < 1 respectively. The 


last inequality follows by applying (4.3). 


In particular, we obtain for 0' = cpR^R 


< ii-q ln-‘I 


A. 


(4.6) 


This type of geometric construction was hrst used by Choksi, Conti, Kohn, and 
one of the authors in |1] in the context of branched patterns in superconductors. 

For the convenience of the reader we reproduce a version of the proof of Lemma 
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Proof of Lemma\4^ We split the proof in two steps. 

Step 1. In the hrst step, we construct a set Qji that covers most of {x = 1} (see 
Claim 1 below) and has radius of curvature < R (see Claim 2, below). As before, let 
Xr = f^R * X denote the mollihcation of x on scale R. We dehne 


^R ■— {Xr > 1/2}- 


This time, we take the non-smooth “Dirac sequence 

f^Rix) = 


^ for 
0 for 


\x\ 

Ixl 


<f 

>f 




so that Qji can be characterized via the density of {x = 1} in balls of radius R/2 as 
follows 

I l{x = 1} n5|(x)| > 1\Br{x)\^ . 

We show now the two following claims. 

Claim 1: for Qr dehned above we have 


[x < R [ |Vx| + [ X- (4.7) 

J J J Qjj 

Indeed, 

[x - [ X = |{x = 1} n {Xij < 1/2}| < 2i|x - XhIIi <2R f |Vx|. 

J J J 

Claim 2: There exists a hnite number N of points Ui G VLr for i = 1, such 

that 

TV 

Qr C [jBR{yi) and ^ ^ ^ j X- 

i=i 2 

Indeed, let {|/i,..., t/tv} be maximal with the property that Bniiii) fl Bniiij) = 0 
for any = 1,...,77, with i ^ j. The hrst part of the claim follows from the 
maximality of {|/i, indeed if there were an yo G Qr with yo ^ BR^yi) and thus 

BR{yo) ni?fi(T/j) = 0 for all i = also the strictly larger set {?/o, ?/i, ?/v} 

would be admissible. 

The second part of the claim can be seen as follows: 

N N 

n’^R-‘ = Yi \Bs{yi)\ < 2^; Kx = 1} n < 2|{x' = i}|, 

i=l i=l 

where in the hrst inequality, we have used that for yi G PIr we have |{x = 1} H 
BR{y)\ > }^\BR{y)\. In the last inequality we have used the pairwise disjointness of 

{B R{yi)}i=i,...,N- 
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Step 2. In the second step, we construct the potential 4>r,l- We introduce the 
capacity potential (j)R^L of -Br(O) in -Bl(O) given by 


iM^) ■= 


( 1 for 

In TV . 



L 

< 

"1 

I 

, L for 

R 

< 

L 

< 



0 for 

L 

< 

L 


J 

1 


£ | 0 , 1 |. 


We dehne 


Claim 3: we have 


i,l{x) ■■= max - Vi) e [0,1]. 

1=1,...,N 


x^R |Vx| + / x(pR,L- 


Indeed in Claim 2 we have seen that VLr C IJ^i BRiVi)-: therefore, since by dehnition 
4>r,l = 1 in 5^(0), we deduce that cjiR^L = 1 on VLr. This, together with Claim 1 


implies Claim 3 and thus (4.2). 

Claim 4- We have 

I 4^r,l<l^r-^ j X- 

Indeed, by the dehnition of 0 r,l, we have 

I <I>R,L <N j ^R,L < NL^ < L^R-^ I X, 

where we have used the upper bound on N established in Claim 2. 

Claim 5: 

J max{-A0/j,L,O} < i?“^(ln"^ J y. 

Indeed, using the well-known fact that the singular part of (—A) max{0i, 02} is neg¬ 
ative, we conclude similarly to the previous step: 


max{-A0R,L,O} < N j max{-A0i{,L, 0} 

C. 


= A27rln-^ 4 < R~^ ( In"^ 4 


R 


R 


This concludes the proof of Lemma 4.2 


□ 


For the convenience of the reader, we reproduce the proof of Lemma 4.1 from [8]. 
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Proof of Lemma 4jJ_' By the three properties of the geometric construction we have 
[ X < R [ | Vx | + [ f)R,LX 


= R I Vxl + / (f>R,L{x - <h) + <h / (fR^L 


< R \Vx\+{ \V<fR,Lf \\v\-\x-n^] +<^> 


l,L 


< 


R / |V\'| 


R-^ I ln-‘ ^ 


\' / l|V|-'(x'-4>)t)’+4(^ 


X- 


We hrst absorb the factor f x of the middle right-hand side term by Young’s inequality 
to get 

fx < /|Vxl + fl-"flii-Yl /l|V|-'(x-<l>)l" + 1>f4') fx- 


In order to absorb the last right-hand side term, we choose L to be a small but order 
one multiple of ^~^R. Since L is a small multiple of we have 1 so 

that indeed we can absorb; since it is an order one multiple of and <h ^ 1, we 

have L 3> -R and In h ~ In Hence we obtain: 

fx < Rj\Vx\ + R-"(in-'P\ J\\V\-Ux-i)f- 

We hnally optimize in R by choosing R = {J |Vx|)“^((ln“^ J ||V 

and we get 


X < ( In ^ ^ 


|Vxl 


ivr'(x-<^>)i‘ 


Dividing by and multiplying by Ins X obtain the desired estimate. 


□ 


As in the previous section, we recall here the weak version of our interpolation 


inequality (1.3) in dimension 2, which was proven in the PhD thesis of Viehmann 


Proposition 4.4 ([IS]). Let u : [0, A]^ —>■ M &e a periodic function with u > —1 and 
fu = 0. 

Then, there exists a constant C > 0 such that 

sup/x(ln3/i) |{|m| > /i}|3 < llVwllf |||V|“^m|||. 

/i>e 
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Proof of Proposition 4^ By Lemma 3A and by a scaling argument in x, it is enough 
to show 

sup/x5(ln5 p) |{|m| > fi}\ < ||Vm||i + lllVr^MlIa- 

fj.^1 

For a given level /i S> 1 we consider the characteristic function G {0,1} of the 

corresponding level set of u, that is 

{Xm = !} = {«> fi}. 

For given length scales R L (to be chosen later) let be the potential con¬ 


structed in Lemma 4.2 based on According to Lemma 4.2 we have 


j Xfi ^ RJ |Vx^| + j 

Using that > 0 and the crucial assumption u > —1, we rewrite this as 

R / |Vx/.l + / x^l(t>^l,R,L 


< R 
= R 

< R 
= R 

< R 

< R 


P P J Xijl4’^i,r,lR 

^Xr.\ + j X|^(|>^^,R,L(u + 1) - j X^^(i>^^,R,L 

'^Xf,\+P~^ / + 1 ) - / Xfi4>fJ.,R,L 


VXul +P 


-1 


,R,lU + P / (1 “ Xti)4’tJ.,R,L 


Vx/^l+h y 4>ti,R,LU + fi ^ j 4>f,,R,L 
^Xl.\ +P~\[ [ ||Vr^M|^ 1 +/i"W (j)^,,R,L■ 


We now insert estimates (4.4) and (4.6) from Lemma 4.2 to obtain 


X^l 


< 


R / |Vx, 


1 

+p-^ ( R~^ ( ln“^ / 11^1”^“!^ 




Xfi- 


In order to absorb the last right-hand side term, we choose L to be a small but order 
one multiple of fi^R. Since L is a small multiple of fi^R, we have P~^(^Y <C 1 so 
that indeed we can absorb; since it is an order one multiple of p,^R and p 3> 1, we 
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have L ^ R and In ^ ~ In fx. Hence we obtain: 


Xfj. 


< 


R I IVXmI + /i M i? ^(In V) / Xm 




In order to absorb the factor f X/i of the last remaining right-hand side term, we use 
Young’s inequality and we get 


|{« >fi}\= X,<R |VxJ +/i-2i?-^(ln~V) 




ivr^M 


(4.8) 


By the coarea formula, we have / \ Vx^l\dxdf^ < J |Vm| so that there exists a 
/i' G with fi f |Vx^/| < 2 J |Vm|. Using (4.8) for jx replaced by /i' we thus have 


\{u > jx'}\ < Rfi ^J\Vu\+ix' ^R ^(In ^/i') 

which because of /i' G [^,/i] turns into 


\{u > fi}\ 


< 




R^^-^ / |Vm|+/ i-"/2-^(ln-V) / llVr^w'^ 


We multiply with /is Ins /i and we get 

4 1 

/iS Ins /i|{M >MI 

< i?(/iln/x)3 / IVm|-|- i?“^(/iln/i)“3 




The choice of i? = (/r In /i) s yields the desired estimate. 

We now give the proof of our strong interpolation inequality in dimension 2. 


□ 


Proof of Proposition 1.2 By an approximation argument, we can assume that u is 
a step function (indeed all the quantities appearing in inequality (1.3) well behave 
under approximation by step functions). By a scaling argument in x and the result 
in Proposition 1.1, it is enough to show for M ^ 1: 


MS Ins u < IIVm||i -h IIIV| 


u\ 


J {u>2M} 

We consider an arbitrary level /i > M 3> 1 and start as in the proof of Proposition 


4.4, considering the potential (t)f_i,R,L- Observe that, since we are assuming that u is 


a step function, then is piecewise constant as a function of /r, and therefore it 

is measurable in /r. This will be important later since we will integrate in d/r. 

/i, we get 


For L chosen such that (4 i" ~ 


w 


Xfj' 


< 


R / |VxmI +/^ / 4 'ii,r,lu. 
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But we now rather proceed as in Proposition 1.1 We multiply with (pln/i)3, choose 
R = (/iln/x)“5 and integrate in p G (M, cx)) for M 3> 1: 


' M 


< 


{fi In /i) 3 / Xfidx d/i 


|Vx^|da;d/i + 


' M 


< liv«lh 


/•“ ln3 /i 

'M p.5 


/•“ ln3 

'M /i5 

4>^l,R,Ldf^ 


4’ti,R,Ldix u dx 


VI ^ullo. 


On the last right-hand side term we argue along the lines of Proposition |1.1 

) of our geometric const 

'^4>ti,R,L ■ '^4>ti',R',udx < 


(4.9) 


now 


using the property (4.5) of our geometric construction, that is 

1 1 


In I 

Jri 


X^dx, 


(4.10) 


where R! and L' are related to jj,' like R and L to fi, that is, R' := (/i'ln/i') 3 , 


(w)^ ~ /i^ Using (4.10) and by the choice of R and L, we get 


In 3 /i 


' M /i3 

poo poo 


4’^i,R,Ldp 


= 2 


l|4.10| 


I I ln3 p ln3 fi' 

M Jm /i3 /i'3 J 

In^ /i ln5 /i' 


Jm Jm p's 

'•00 rfl ^ ^ ^ 


< 


'M Jm /i3 p's 

r r [ x.dxdf, 

'm Jm /i'3 J 

poo p 

/ (/ilnp)5 / Xfidxdfx. 

'M J 


'^4>fi,R,L ■ '^4>^i',R',L’dxdp'dp 
■ V(j)^/^R>^L'dxdfi'dfi 
Xiidxdfi dfi 


Hence, coming back to (4.9), we can absorb this term by Young’s inequality and 
obtain 

/ (/ilnp)^ / XM^2;d/i < ||Vm||i + |||V|~^n|| 2 . 


IM 


We conclude by observing that for M 3> 1 


IM 


(/ilnp)3 / Xfj.dxdiJ,= 


r'u(x) 


i{u>M} Jm 
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(/i lnp)3d/i dx> 


^ 1 - 
M3 Ins u. 


'{u>2M} 

























□ 


5. Proof of Propositions 11.31 and 12.11 
In this section we prove Propositions 1.3||2.1 , and 1.4 The two main ingredients in 


the proofs are the geometric construction of Lemma 4^ and the Kantorovich duality 
for the Wasserstein distance. 


Remark 5.1. In the proof of Proposition 1.3 we need the analog of the geometric 


construction of Lemma 4.2| in any dimension d. Following the proof of Lemma 4^ it 
is easy to see that given a function y : [0, A]'^ —)■ {0,1} there exists a set VLr, which is 
dehned as 

I I{X = 1} n5|(x)| > \\Br{x) 

and N points yi G VLr for i = 1,..., such that 


N 


C U BR{yi) and ^ 


A, 


(5.1) 


2 = 1 


where {i/i,..., y^} is maximal with the property that 5i?/2(i/i) ni?R/ 2 (i/j) = 0 for every 
i,j = 1,..., N such that i ^ j. 


Similarly to Lemma 4.2, we want to dehne now a potential associated to y. 
Since here we are not interested in the logarithmic behaviour of the potential, it is 
enough to dehne (j)R as the characteristic function of With this choice 


of (l)R, similarly to (4.2), we have 


X $ / |Vxl + / X>I>R- 


(5.2) 


We are now ready to give the proof of our Proposition 1.3 


Proof of Proposition \1.3\ As before, by approximation we may assume that m is a 
step function. By a scaling argument in x it is enough to show that there exists a 
constant C only depending on d such that 


ll(« - < C (IIVulli + Wl(n,l)) 

3d 

To make the proof more readable, we divide it in three steps. 


(5.3) 


Step 1. We start as in the proof of Proposition |4.4| using our geometric con¬ 
struction. For a given level /i, let as before ^ {0; 1} denote the characteristic 

function of the set {u > y} and let be the potential associated to y^ dehned in 
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Remark 5.1 By (5.2) we have 


Xfj. 


< 


R / \^Xf,\ + / 


R 


- R J R ~ J 


(5.4) 


Let e be a small parameter to be adjusted later. We multiply (5.4) by p(2+3d)/(3d)^ 


we choose R = C\Rjj, ^/(srf) (-^v^here C\ is a dimensional constant to be specihed later) 
and we integrate in J — for /i > 1/e'^, to get 


r' + OO 


+ 00 


L/e'* J h Jl/e<^ 


|Vx^|(ixd/i 


'•+00 


dfi 


/ — ]u{x)dx. 

n/sd- h 


(5.5) 


Using the coarea formula as before, the hrst term on the right-hand side is estimated 
as follows 


b+OO 


h/ed 


\Vxp\d^d^ ^ l|Vn||i. 


To estimate the second term on the right-hand side of (5.5), we set 

/• + 00 


^p{x) := / ^i3d(j)^^^(x) — . 

Jlle^ 


(5.6) 


Using the Kantorovich duality (1.5) with n = 1, with replaced by —'ll!, and with 
cost function c(x) = |x — we have that 


^p{x)u{x)dx < -^W^iu, 1) + / ^{y)dy, 


where 


\x — y\^ 

%lj{y) := sup <! (p(x)-r— ^ = sup 


r^+oo 


h/ed 


y^'^(Pti,R[x 


dfi \x — I/p 

- 9 - 

/i 


(5.7) 


(5.8) 


Combining all together in (5.5) we get 
h 


' 3 ^"" [ Xi^dx— < ||Vm||i + ^W^{u,l) + [ 'ip{y)dy, 


n/ed 




(5.9) 


where ijj is dehned in (5.8). 


Step 2. In this step, we estimate the term J 'ijj. We will show that, for £ small 
enough, it can be absorbed on the left-hand side and this will conclude the proof. In 
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this step, we will assume that the following inequality holds: 


< sup 

X 

<2 I 


2 dfi \x — y\‘^ 

'Ijs^ /X S 


+00 


' sup <j 

l/e<i 


\x — y\^ I dy f~^°° ~ dfi 
— =■ 

4- h Jxje'^ y 


(5.10) 


where, for simplicity of notations, we write sup{/}+ in place of (sup{/})_,_. We will 

. B 

dy 


prove this inequality in Step 3 below. By (5.10) we have that 

f*+oo 


i^{y)dy < 


/ i^^l{.y)—]dy = 

'Ije^ /X 


^+00 


' l/e^ 


dy 


i^f2{y)dy ] 

) y 


(5.11) 


We recall that (see Remark 5.1) for any /i, is the characteristic function of the 
union of N balls Bn^yi) for i = I, where N is bounded above by ^ 

This implies that 

k - y\^ 


ij^iy) = 2snp<Ciy'-id(j)^ji{x) 


= 2 


2 

Cifi^ if yEBn^yi) for some i = 1,..., iV 


max 


jCi/iSd - dl/e'^, o| ify ^ BR{yi) for every i. 


where dy := dist (^y, U^i BR^yi)^ , with dist{y, A) denoting the distance between y and 

the set A. Observe that if y is such that dy > Ci then ^^^.{y) = 0. Moreover by 

the choice R = C \^‘^and by our assumption y > l/e'^, we have > R. 

Thus 'ifjy is supported in the union of N balls Bi{yi), with radius I = R + ey^ < 

2ClBsyJd. 

Hence we have 


i^t2{y)dy < Ny^Rd 


^ I Xt^dx < eV 


2+3tf 


X^^dx. 


Using (5.11) this yields 

/ i^{y)dy < 


A fed 


^f2iy)dy ) y ~ 


l/ed 


y 


2+3d 

3d 


dy 

Xfddx — . 


Combining all together in (5.9) we get that there exists a constant C > 0 only 
depending on d such that 


'l/e' 


y"^ I X,dx^ < C (^11 Vnlli + ^lU|(n, 1) + 
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f°^ 2+3d f dy 

/ y^d Xf,dx — 

Ayd J y 









Choosing e = (2(7) we can absorb the last term on the right-hand side to get 


2+3d 




I2C 


J Xfidx— < 2C 


1 + {2Cf/'^W^{u, 1) 


Evaluating the integral in /r on the left-hand side we deduce (5.3) with C = 4C^. 

Step 3. In this last step we show (5.10) and therefore we conclude the proof. By 
dilation in /i, it is enough to prove that there exists a dimensional constant Ci such 
that 


sup 


^+oo 


2 _ 


dfi \x-y\‘^ 

-^— 

fi 


b+oo 


< 2 


+ 


sup - 


\x — 1 d/i 


/i 


(5.12) 

We start by proving two technical estimates. 

Claim A. Given a function /, we have the following relations between (discrete) dyadic 
sums and (continuum) logarithmic integrals: 


1 

2 


„2 r+oo j j‘2 °° 

I ^/(«2‘)<i»< / f(i^C< / Y.!{e2’‘)de. 

1 fe = 0 h- Jl 


(5.13) 


The proof of Claim A is trivial; indeed, using the change of variable fi = 62^, we have 

< 2 


r-2 °° 


oo 

'£,/(D2'’)M = J2 I, f(D2’‘)M = Yl 

k=0 


k=0 


k=0 




A 


> 


/i 

d/i 
/ h — 
h 


Claim B. For any fixed d G M and k G N, let (pe^k be a characteristic function. 
Then, for p > 0, we have that for any x the following estimate holds: 


2P 


^{e2AA,kA < sup {{e2AA,kA}- 

^ ^ 0<fc<oo 


k=0 


To prove this claim, we set 


K{x) := sup{A:| AA^) 7^ 0}- 


(5.14) 


(5.15) 


If K{x) = oo inequality (5.14) holds trivially since the right-hand side is inhnite. 
If K{x) < oo then we have 

J2i02AAAA < J2(e2A = o^ 

k=0 k=0 

,p 2P 


2P-1 

2P 


£ ^A-r = sAt sup {(92‘)!>.#.ii.»(n}. 

ZP — i ZP — I 0<k<oo 


which concludes the proof of Claim B. 
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Now, using (5.13) and (5.14) we deduce that (5.12) holds. Indeed, recalling that 
is the characteristic function of a hnite union of balls (depending on /r), we have 



where Ci = y** (cf. Claim B). 

233-1 


This concludes the proof of (5.10) and thus of the proposition. 


□ 


We give now the proof of the interpolation estimate in additive form stated in 


Proposition 2.1 


Proof of Proposition 2^. We follow the same strategy of the proof of Proposition |1.3| 
and we divide again the proof in three steps. 


Step 1. We start by showing that inequality rt2.21| holds with = 1, that is we 


want to prove that for every u, n > 0 such that A J u = A = there exists 
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a constant C with <h < C* snch that 

/ 3d+3 / X 

{u - l)r^ < C (||Vnlli + Wi{u,v) + \\\\/\-Hv - $)i|2j . (5.16) 


We begin as in the proof of Proposition 1.3, assnming that n is a step fnnction and 


nsing onr geometric constrnction. We have (cf. (5.4)) 

j i^wi + 


(5.17) 


3d+3 


. -. I u / ]^/2 2 

We mnltiply (5.17) by /i^d+i, choose R = C 2 /i ^d+i (where C 2 is a dimensional 
constant to be specified later) and integrate in J ^ for /i > 1, to get 


3d+3 


/ Xix{x)dx— 
jl 


< 


\V x^l{x)\dxd|Ji + 


) u{x)dx. (5.18) 

li 


Using the coarea formnla as before, the hrst term on the right-hand side is estimated 
as ^ 

j\'^X^^{x)\dxd|x<\\Vu\\l. 

To estimate the second term on the right-hand side we proceed as in the proof of 
Proposition 1.3 By analogy with ( |5.6 ), we set 

djx 


ip{x) := / /r 


2 

3d-|-l i 






Using again the Kantorovich dnality, we obtain 

/ ip{x)u{x)dx < !I) + / i’(y)v(y)dy. 


(5.19) 


where analogously to (5.10) we have 

= sup{(p(x) - |x - I/p} 


< sup 


^+00 


djx 


- \x-yY 

/i 


'‘+00 


< 2 


s^^i\C2^^^'^+^(|)^x,R{x)-\x-y\ \ — 


dfj, 

y 


where Co = 


2 

2 3d+l 

2 • 

23^-1 


Plugging (5.19) into (5.18), we get 


3d+3 


dfi 


y Xix{x)dx—<\\Vu\\i + W 2 {u,v) + J ij{y)v{y)dy. 
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(5.20) 













Step 2. In this step we estimate the term f 'ipv. In order to do that, we define the 
functions 

and ^{y) := [ —• 

Ji h- 


ijf.iy) := 2sup C'2/U3d+i0^^^(x) - \x - y\^ 


The second term on the right-hand side of (5.19) is bounded by 


ipv < J ^pv = J 'ip{v — $)-!-$ j ip 

< II |V|j?||2|||V|-^l.-<*)||2 + ® f 


(5.21) 


We give now an estimate for the quantity || |V| 2 'i/'l| 2 - We recall that here refers 
to the homogeneous fractional Sobolev space, endowed with the seminorm dehned by 

\f{x)-f{x)\^ 


|V|5/||2(M'^) : = 



I cc cc I ^ 


-dxdx. 


We have 




r 7 dy 


'1 



fOO /■oo 


'1 Jl 



2 rdy rdfi 

Jl Jl 


^ rdy r 

Jl ii /^' 


dydy'\ 

\y-W^^ y y' ) 

{^^{y) - ^f,{y)){^^>{y) - t/V'(^)) 


\y-y\ 


d+l 


dydy 

dydy 


where in the last equality we have changed the order of integration and we have 
used the symmetry between p and y'. To estimate the quantity I^'^' we proceed as 
follows. First recall that, since is the characteristic function of the union of N 
balls Bji{yi), we have 


i’l^iy) = 2 [snp<C2iJ,3‘^+^(j)i,,R{x) - \x - yp 


= 2 


C2/i3<i+i if ?/ G Bji{yi) for some i 

max I(72/13d+i — 0 I iiy^Bfi{yi) for every h 


(5.22) 


where dy := dist{y, IJ^i Observe that if y is such that dy > then 

ipniy) = 0, and that by the choice R = we have y^d+i > for y > 1. 
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^^ 1 /2 ^ 
Thus is supported in the union of N balls Bi{yi), with radius I = R + C 2 < 

1 /2 1 

2 C 2 /i^a+T. Moreover recall that the number N of balls is bounded by 


^ ~ p^d I 


By (5.22), we see that ipa satishes 




2 

^ n3d+l 1 


(5.23) 


(5.24) 

(5.25) 


Using the fact that -0/^ is supported in LiiLiBi{yi) and the symmetry between y and 
y, we have 


< 2 V / 

i=l 
N ^ 

= 25: 




dydy 


=1 dBiivi) J{y.\y-y\<r} 
N , 


i^f^jy) - ^yiy))i^^.'iy) - ^y'jy)) 

\y-y\d+l 

i^^.iy) - ^f.{y)){^y'{y) - 


i=l d Bi{y,) J {y:\y-y\>r} 


\y-y\ 


d+l 


dydy 
dydy 


N 




2=1 


To bound the hrst term , we use that I < gradient bound (5.25) for 

and spherical coordinates centred at y. We get 

f f i^yjy) - (y)){^y'iy)- (y )). 

'Bi(yi) J{y-.\y-y\<r} 


Ty,y _ 




< 


\Biiyi)\ llVV'^.lloo W^'ip, 


\\ 00 


\y-y\ 

r 

L 


-dydy 


dp 


< 


^ 3d+l jP 3d+l f = jj 3d+l ^ 3d+l 3d+l 


On the other hand, using now the L°°-bound (5.24) for ipu, we deduce 


Ty,y _ 
-'*,2 — 


< 


{'4’y{y) - '<PM){'>Py'{y) - ^u(2/)) 


JBi(yi) J{y.\y-y\>r} 

IIV^Atlloo llV^/.i'lloo 




\y-y\ 


dydy 


r-oo pd-l 


P' 


d+1 


dp 


< rp3d+l Id 3d+l _ = //3d+l //3d+l II 3d+l 


p3d+i i^ad+i n 
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Optimizing the sum in r , we get for r = p2{3d+i) ^ 2 {id+\) 

N 

< 2^(/a'‘' + /5'‘') <]V^! 


2=1 

1 2rf+3 ^^ 


2d+3 , 3 

2(3d+l) jj 2(3d+l) 


X f o I 2d 2d+3 , o i 

< —/i2{3d+l)^ 2{3d+l) / < p3d+1^2(3d+l)^ 2{3d+l) / 


where we have used the bound (5.23) for iV and the choice R = 2/(3rf+i)_ 

Finally, integrating in /i and jj,' we get 

^ Ji d-' 

f°° dfjL dfx' 6d+3 ,_3 _ f 

< I Jl i _lF^2(3d+i)^ 2C3d+i) dxx^l (5.26) 


| V | 2 F ;||2 = 2 / 


< 


'1 d' JI 

f°^ 3d+3 f dp 

I y 


We now estimate the second term $ / V' appearing in (5.21). By an analogue but 
simpler computation, we deduce that 


"" / N ^ ^ A ^ —A ^ A I z a i 3d+2 # 

i^t.{.y)dy < N< R Vsd+l/'i I _ ^3d+1^3d+l;^3d+l / ;^^ = ^3d+l / 


Thus, plugging this into ^ f we get 




pirr J X^ddx^. (5.27) 


Plugging (5.21), (5.26), and (5.27) into (5.20), we deduce that there exists a constant 


O > 0 such that 

^“pi^ J ^^^<C\\Vu\U + CW^{u,v)+ 


(5.28) 


3d+3 / , ap \ 3d+2 f , dp 


+ 0|^|||V| 2(u-<F)||2 J J y^^dx^^ y 

(5.29) 

Since (3d + 2)/(3d+1) < (3d + 3)/(3d +1), we can absorbe the last term on the right- 
hand side for d) < ^ . Moreover, using Young inequality, we can absorb the second 
factor of the third term on the right-hand side. Finally we evaluate the integral in p 
on the left-hand side to get 

/ 3d+3 / N 

(u - l)r^ < 20 (llVulli + lF|(u,u) + WM-Hv - $)||2j . 
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(5.30) 


















Step 3. In this last step we show that, by scaling, the following inequality in 

3d+l 

additive form containing the parameter v holds for d) < z/3‘^+3 /{^C) 


3d+l 
U — Z/3d+3 


3d+3 

3d+l 


< 2C [ IIVmII- 






Indeed, set u = Mu, v = Mv, d) = M<h, and x = Lx, by (5.30) we have for 
$ < (2C'M)-i 




3d+3 


. 3d+3 


We divide by L“M 3^+1 to get 

3d-|~3 ^ 2 


y oa-|-o 2 / ^ 2 3d_1 1 ^ 

<L-^M-33+t / |VM|+L2M-33+Tpp|(w,{))+LM33+T|||V|-2(h-<h)|| 
2 

Choosing L = M 3d+i we have 


3d+3 


(u — M < / IVmI+M 3 <i+ipp'|(.y^{))-|-]\^3d+i |||v| 2 ({) — <|))||2_ 


J-1 3d+3 

This implies ( |5.31[ ) for i/ = M 3d+i. □ 

Finally, we prove Proposition [L^ which follows easily by Proposition |2.1| by scaling 
arguments. 


Proof of Proposition l.f As a corollary of Proposition 2.1, we deduce that given two 

bd 


nonnegative functions u, v : 
additive form holds 


-)■ 


with f u = J V < 00 , the following inequality in 


3d+l 

U — i/3d+3 


3d+3 

< ||VM||i + z/3^hF|(M,t;) + z/^|||V|-W||2. 


We pass from the additive form to the multiplicative one, by scaling in x. Indeed for 
X = Lx, we have 


3d+l 

U — Z/3d+3 


<^-Nlv«lh 


2d 

+ Ld+i 


vL) 'i+i IF|(m, v) + (z/L) d+i III V| 


1 

'^v\ 


Setting u = uL and choosing L 
we deduce 


(||V«||.)’"‘ 



d+1 

3d+l 


3d+l 
U — i/3d+3 


3d+3 

3d+l 

+ 


;= (llVfiii) 


2d 

3d+l 




d+1 

3d+l 


Taking the supremum in v (that is the supremum in z/ on the right-hand side, since 
z/ = Lv) and raising to the power (3d + l)/(3d + 3) we conclude the proof. □ 
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